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Memory-induced anomalous dynamics: emergence of diffusion, subdiffusion, and 
superdiffusion from a single random walk model 

Niraj Kumar, 1 Upendra Harbola, 1 and Katja Lindenberg 1 

department of Chemistry and Biochemistry, and BioCircuits Institute, 
University of California San Diego, La Jolla, California 92093-0340, USA 

We present a random walk model that exhibits asymptotic subdiffusive, diffusive, and superdiffu- 
sive behavior in different parameter regimes. This appears to be the first instance of a single random 
walk model leading to all three forms of behavior by simply changing parameter values. Further- 
more, the model offers the great advantage of analytic tractability. Our model is non-Markovian in 
that the next jump of the walker is (probabilistically) determined by the history of past jumps. It 
also has elements of intermittency in that one possibility at each step is that the walker does not 
move at all. This rich encompassing scenario arising from a single model provides useful insights 
into the source of different types of asymptotic behavior. 



PACS numbers: 02.50.Ey,05.10.Gg,05.40.Fb 

I. INTRODUCTION 

The use of random walk models in statistical physics 
dates back to the very earliest days of the subject [lHa]- 
They have been used so broadly and pervasively that it is 
impossible to make proper reference to so big a subject. 
The original random walk models were typically ones in 
which the walker moves via a series of random transitions 
characterized by finite length scales for each step and fi- 
nite time scales between transitions. As long as these 
general features characterize the walk, the mean square 
displacement of a walker from its point of origin in the 
absence of a bias, calculated as an average over many ran- 
domly generated trajectories, grows linearly with time. 
If there is a bias, then the mean square displacement 
around the average trajectory grows linearly with time. 
This linear growth has become the universal identifier of 
what is known as "normal transport." Diffusion is the 
quintessential macroscopic normal transport mechanism, 
and is often arrived at by taking appropriate long-time 
and short-distance limits of random walks [J, [7| • 

It is then not surprising that any random walk (or, 
for that matter, any stochastic process) that leads to a 
mean square displacement that does not grow linearly 
with time is called "anomalous." This characterization is 
shared by stochastic processes whose mean square dis- 
placement grows either sublinearly or superlinearly with 
time, the former often called "subdiffusive" and the lat- 
ter "superdiffusive" processes. In the world of random 
walks, subdiffusive behavior is often associated with wait- 
ing time distributions between steps that have fat tails, 
e.g ., that decay as an inverse small power of time [8|- 
Il0| . The average time between transitions in this case 
is infinite. On the other hand, superdiffusive behavior is 
typically associated with step length distributions that 
have fat tails, so that the average distance per step is in- 
finite [lfj| | . Levy flights are well-known examples of such 
processes [10| . 

In recent years there has been considerable interest in 
formulating stochastic models that can exhibit different 



types of behavior in different parameter regimes. More 
specifically, for example in the context of random walks, 
it has been observed in a variety of contexts that the 
behavior of a walker may be normal in some parameter 
regimes but anomalous in others. This is almost trivial 
to envision in terms of the quantities already introduced 
above. For instance, consider a model with a waiting 
time distribution with a power law tail of the form t~ a . 
It may then happen that in some regimes of the model the 
waiting time distribution decays slowly (a < 2), so that 
the average waiting time between transitions diverges. In 
a different regime of the model the decay of the waiting 
time may be sufficiently rapid (a > 2) for there to be a 
finite mean transition time. The process would be subd- 
iffusive in the first regime and normal in the second. On 
the other hand, for example if a model has a jumping 
pattern with an infinite average jump distance in some 
regime but a finite one in another, the model would ex- 
hibit superdiffusive behavior in the first and normal dif- 
fusive behavior in the second. Transitions between nor- 
mal and anomalous behavior within a model can occur in 
many other ways as well. The literature is far too large 
to reference it properly. Perhaps the earliest of these 
in the random walk context was introduced by Zumofcn 
and Klaftcr [11] to describe the dynamics generated by 
iterated maps. A more recent example directly related 
to the discussion in our paper deals with a class of ran- 
dom walk models called "elephant random walks," a term 
coined by its creators [12j because they involve a sort of 
perfect memory often (and most likely inappropriately) 
associated with elephants. In these models, which have 
also been mentioned in applications as diverse as ecology, 
economic data, and DNA strings, the memory induces a 
transition between normal and superdiffusive behavior 
observed with a change in parameter values |13l4l9| . 

While there is thus a plenitude of stochastic models in 
which the mean square displacement changes from one 
behavior to another as a parameter of the model is var- 
ied, it is more difficult to find models that exhibit all three 
forms of behavior, that is, ones in which a change of pa- 



rameters leads to asymptotic subdiffusive, normal, and 
superdiffusive behavior. Such versatility can be found 
in generalized Langevin equations [20J and in dynam- 
ics governed by fractional Brownian motion or fractional 
Langevin equations (a helpful list of references to the 
fractional Langevin literature can be found in 21]). How- 
ever, to the best of our knowledge, there is no single 
random walk model that exhibits all three forms of be- 
havior. In this paper we present such a model. It is a 
random walk model with a memory, and a simple sweep 
in parameter values indicating the direction of motion of 
the next change can lead to transitions covering all three 
forms of behavior. The model is inspired by the "ele- 
phant random walk" introduced in [121 ] . In Sec. ITTl we 
present the random walk model. In Sec. lIIII we calculate 
the moments of interest to characterize the nature of the 
motion of the random walker. A detailed analysis of the 
results is presented in Sec. IIVI We conclude with a brief 
recap in Sec. [V] 



II. THE MODEL 

Following the notation in [12[, consider a random 
walker on a one-dimensional infinite lattice with unit dis- 
tance between adjacent lattice sites. Steps occur at dis- 
crete time intervals. At each step the walker can take 
one of three actions: it can move to the nearest neighbor 
site to its right, to the site on its left, or it can remain 
at its present location (that is, the walk is intermittent). 
We denote the position of the walker at time step t as Xt- 
The position of the walker at time step t + 1 is 



X t+ l = X t + (Jt+l- 



(1) 



Here <7t+i is a random number which can take on one 
of the values —1, or +1. The choice of this random 
number at each step depends on the entire history of the 
walk, {<7 t } — (ci, ...,<Jt), as follows. A random previous 
time k between 1 and t is chosen with uniform probability. 
If o~k = ±1, with probability p the walker takes the same 
step at time t + 1, i.e., <7t+i = <7fc. With probability q, 
the walker takes the opposite action, o~t+i = —o~k- The 
walker can stay at rest with probability r, and of course 
p + q + r = 1. If o"fe = 0, the walker stays at rest with 
probability 1. The process is started at time t = 1 by 
allowing the walker to move to the right with probability 
s and to the left with probability 1 — s, i.e., the first step 
excludes the possibility that the walker may not move. If 
the walker is initially at X = 0, then the position of the 
walker at time t > 1 is given by 



Xt 



E 
fc=i 



Cfc. 



(2) 



We make a special point here about the significance of the 
parameters p and q: they are not the familiar parameters 
that indicate a next-step asymmetry. Rather, they are 



memory parameters that indicate whether the walker, if 
it moves at all at the next step, is likely to follow the 
randomly chosen past step (persistence probability p) or- 
is, instead, a rebellious walker who does the opposite 
(probability q). It is thus not straightforward to a priori 
predict the direction and growth properties of the walk 
as a function of these parameters. For r — 0, the model 
reduces to that of (12| . However, we will show that the 
inclusion of a probability that the walker neither follows 
nor rebels against the randomly chosen prior step, r ^ 0, 
turns out to be the crucial generalization that leads to 
the all-encompassing model. 



III. MOMENTS 

The quantities of interest to characterize the nature 
of the motion of the walker are the mean displacement 
and the mean square displacement as a function of time, 
both of which we are able to calculate analytically. The 
parameters of the problem are two of the stepping proba- 
bility parameters; those that turn out to be most useful in 
characterizing the properties of the walk are the "staying 
probability," r, and the memory asymmetry parameter, 



7 =P~Q- 



(3) 



We start by noting that for a given history {<r t }, the 
conditional probability that ot+\ — c, where t > 1, can 
be written as, 

P[<j t+1 = <j\{a t }} = 1 - a 2 



2t 



fc=i 



and for t = 0, 



P[°i =<r] = -(l + (2s - 1)*) 



(5) 



Using Eq. (j4}, the conditional mean values of at+i for 
t > 1 in a given realization is given as 



(o-t+ilWt}) = Yl aP \- a t+i = a \{ a t}] = -7 x t, (6) 

cr=±l,0 



which, on averaging over all the histories, gives the fol- 
lowing mean values, 



\<rt+i) 



7 



xt), 



(7) 



where (xt) is the mean displacement of the walker. Per- 
forming the average of Eq. (QJ and using Eq. ([7]) leads to 
the recursive equation 



(xt+i) 



1 + £)<•<■/}. 



(8) 



whose solution is Ili 



(x t ) = (2s -1) 



r(i + 7 )r(t) 



f for i > 1. 



(9) 



The mean position of the walker vanishes if s = 1/2. 
For s > 1/2 and s < 1/2, the mean position is positive 
and negative respectively. Thus, the first step, and the 
first step alone, determines whether the walker moves to 
the right or left macroscopically. The further evolution 
of the mean displacement does depend on the parameter 
7. We note that for a symmetric memory, 7 = (p = 
q), the mean position is time independent. For 7 > 
(p > q), the mean position increases with time with an 
exponent which is smaller than unity for nonzero values 
of the rebellion parameter q, and so the velocity of the 
walker decreases with time. For 7 < (p < q), the mean 
position decreases monotonically with time at long times 
and thus, on average, the walker returns to the origin 
and thus remains localized in space. This type of motion 
might be useful to model the dynamics of animal home 
range behavior [22j. 

We next compute the second moment of the displace- 
ment, (x 2 ). For this, using Eq. ([T]), we note that 



xLi = x 2 + a?j_i + 2x t a t+ i. 



u t+i 



't+i 



(10) 



For a given history {at}, the conditional average of x 2 +1 
is then 

<a£ +1 |{<7 t }) = x 2 t + (a 2 t+1 \{a t }) + 2x t (a t+1 \{a t }). (11) 

Here we have made use of the fact that for a given 
history {at}, Xt is known, which allows us to write 
(x t a t+ i\{a t }) = x t (a t+1 \{a t }). Using Eq. g]), the condi- 
tional mean values for a 2 +1 can be written as 



which on averaging over all histories gives the recursive 
relation 



(a 2 +1 ) = (l-- t )(a 2 ), 



leading to the solution 



2v _ r(t-r) 



r(i-r)r(t) 

Using Eq. HT5J) in flU]) gives 



\k=i 



T{t-r + l) 
r(2 - r)T(t) ■ 



(15) 



(16) 



(17) 



Substitution into Eq. (I13[) and solution of the recursion 
equation leads to the mean square displacement for 7 / 



<*?) 



1 



(2 7 + r _i)r(t) V r(2 7 ) 



r(i + 2 7 ) r(i + t-r) 



r(i - r) 



1 



( 



f2- 



(2 7 + r- 1) \ r ( 2 7) r(l -r)J ' 



(18) 



Note that the dominance of the first or second terms is 
separated by the line r — 1 — 27. For 7 = 0, we find 



<*?> 



1 



r(i + t-r) 

T(2-r)r(t) ~ r(2-r) 



f 1 



(19) 



Equations ©, (TT5]) and (|T^|) are the central results of 
this section. For r = 0, these equations reduce to those 
obtained in [12| • We now proceed to analyze the various 
regimes of behavior of the variance 



Var = (x 2 t ) - (x t f 



(20) 



embodied in these results when r ^ 0. A phase diagram 
associated with this discussion is shown in FigQ] 



fc=i 



(12) 



Using Eqs. (|7J| and (IT21 in Eq. pT|) and averaging over 
all histories, we arrive at the recursive relation 



2 J 
t 



(xl +1 )=[l + ^)(x 2 t)+ J ' 7 ^~ 2 




(13) 



For r = 0, a\ is always 1 and so the term (J2k=i a k) ^ s 
simply t. For r = 1 the sum is unity because only the 
first step contributes. However, for nonzero values of r, 
a 2 , can be either 1 or 0. 

In order to compute (X)fe=i cr fe); using Eq. (fT2|) we 
write 

(*?+ilW) = ^Vt l{^-i}l> + ^^l (14) 



IV. ANALYSIS 

The most interesting result of this calculation is, of 
course, the occurrence of all regimes of behavior, subdif- 
fusive, diffusive, and superdiffusive. However, there are 
additional points to be emphasized because they are not 
necessarily intuitively obvious. 

1. For 7 = or p — q (unbiased memory), the mean 
square displacement increases sublinearly with time since 
the exponent 1 — r is always less than unity for nonzero 
r values. In this case, the mean displacement is indepen- 
dent of time, and so the variance increases sublinearly. 
The behavior is therefore subdiffusive except at r = 0, 
where it is diffusive. This is an interesting, perhaps even 
counterintuitive result, which says that even the small- 
est probability of remaining at a given site at each step, 
with no other ostensible asymmetry toward the site first 
stepped on, leads to subdiffusive behavior. 




FIG. 1: (Color online) Phase diagram for the long-time be- 
havior of the variance Var of the random walker. The left- 
most triangle, labeled as Sub(i 1_r ) (light blue in color), indi- 
cates one regime of subdiffusive behavior. The next triangular 
wedge, Sub(i 27 ) (yellow in color) is also a regime of subdif- 
fusive behavior because here 7 < 0.5 so the time exponent is 
less than unity. The two are separated by the line r = 1 — 27 
(p — 3q). The last triangular wedge, Sup(t 27 ) (light green in 
color) is a superdiffusive regime because here 7 > 0.5. The 
solid line connecting the subdiffusive and superdiffusive re- 
gions (red in color), as well as the r — line for 7 < 1/2 (also 
red in color), indicate the regimes of normal diffusive behav- 



ior because here the variance grows as Var 



+ 2", 



t. The 



inclined solid and dashed thick lines indicate the variation of 
the behavior of the variance for fixed values of the persistence 
parameter p. The dashed (blue in color) line is for p = 0.3 
and, as q and r vary for this fixed value of p, the walker's 
behavior remains subdiffusive throughout, although the form 
of the exponent changes. The solid (also blue in color) line 
is for p = 0.625. Now as q and r vary, the walker's behav- 
ior covers all three regimes of behavior: subdiffusive (of two 
sorts), diffusive, and superdiffusive. 



2. For any fixed 7 < we noted earlier that the mean 
position, which is necessarily one site away from the ori- 
gin at the first step, decreases with time. The mean 
square displacement increases sublinearly with time as 
t 1 ~ r when r > 0. The variance of a rebellious walker 
who has even the smallest probability of staying at a site 
at each step is thus simply subdiffusive. 

3. For fixed 7 > the behavior is more varied. For 
< 7 < 1/2 the mean square displacement grows sublin- 
early with time, as does the square of the mean displace- 
ment (if it is not zero to begin with) . The resulting vari- 
ance thus grows subdiffusively for all values of the param- 
eters subject to this condition, even though the walker 
more often than not follows its previous randomly chosen 
step. However, there are two distinct forms of subdiffu- 
sive growth, separated by the line r = 1 — 27 (p = 3g). 
These two distinct forms correspond to the dominance of 
the first or second terms in Eq. (fT5)) . In one, subdiffusion 
is again caused by the possibility at each step that the 
walker may not move. In the other, the walker's rebellion 
pulls it back. 

4. For 7 = 1/2 and r ^ the growth of the mean 



square displacement is dominated by the first term in 
Eq. (|T5J|, ~ t. For 7 < 1/2 and r = it is dominated by 
the second term, and again grows as ~ t. Together with 
the contribution of the mean displacement, if present, 
this leads to a variance that grows linearly with time, 
that is, the motion is diffusive. The diffusion coefficient 
is given by 



D 



1 



2(2s-l) 2 



3-4p 



(21) 



5. The point 7 = 1/2, r — (where the two red lines 
meet in color) the motion is marginally superdiffusive, 
that is, the variance grows as Var ~ tint. 

6. When 7 > 1/2 the behavior is superdiffusive. Both 
the mean displacement (if s 7^ 1/2) and the mean square 
displacement contribute to this behavior of the variance. 

7. The analysis is considerably more interesting if in- 
stead of following constant 7 lines, as we have done above, 
we follow the behavior at constant r, allowing p and q 
to vary. For small r the figure shows that we cover all 
regimes of behavior as p increases relative to q (7 in- 
creased), namely, the two forms of subdiffusive behavior, 
diffusive behavior, and superdiffusive behavior. However, 
if r exceeds the value 1/2, then we observe only subdif- 
fusive behavior - the tendency not to move is too strong 
to allow anything else. 

8. The analysis is also very interesting if instead we 
follow the behavior at a given value of p, letting q and r 
vary. This is shown by the two inclined lines in the figure. 
The dashed line is for p — 0.3 and the behavior remains 
subdiffusive throughout. The persistence is too low to 
allow for superdiffusive or even diffusive behavior. On 
the other hand if the persistence parameter is sufficiently 
large, as in the solid inclined line (p = 0.625), we sweep 
all of the behaviors as q and r are varied subject to this 
constraint. The value that separates the two regimes is 
p = 1/2. 



V. CONCLUSION 

We have thus presented a random walk model with 
memory that exhibits all three forms of asymptotic be- 
havior as the parameters of the model are varied. The 
characterization has been analytical, in terms of the first 
and second moments of the motion. The model, which is 
inspired by one introduced earlier [ljl [l3| j has only three 
simple parameters. One is a parameter that character- 
izes the very first step of the walk (s). A second is the 
probability that the next step of the walk copies the di- 
rection of a randomly chosen earlier step (p). The third is 
the probability that the next step performs the opposite 
motion as that of the randomly chosen earlier step (q). 
A fourth, which is actually the crucial new parameter of 
the model, is the probability r that the walker simply 
does not move at the next step. It is constrained by the 
others via the conservation condition p + q + r = 1. 



The model can be extended in a number of directions. 
For instance, one can explore changes in the nature of 
the memory. It is also interesting to calculate quantities 
other than the first two moments, perhaps even the full 
distribution and statistics of extrema. Perhaps more in- 
teresting would be an adjustment of the model so that 
it can provide insights into real behaviors that exhibit 
these different regimes. 
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